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Abstract 

Species introductions to new habitats can cause a decline in the population size of competing native 
species and consequently also in their genetic diversity. We are interested in why these adverse 
effects are weak in some cases whereas in others the native species declines to the point of extinction. 
While the introduction rate and the growth rate of the introduced species in the new environment 
clearly have a positive relationship with invasion success and impact, the influence of competition 
is poorly understood. Here, we investigate how the intensity of interspecific competition influences 
the persistence time of a native species in the face of repeated and ongoing introductions of the 
nonnative species. We analyze two stochastic models: a model for the population dynamics of both 
species and a model that additionally includes the population genetics of the native species at a 
locus involved in its adaptation to a changing environment. Counterintuitively, both models predict 
that the persistence time of the native species is lowest for an intermediate intensity of competition. 
This phenomenon results from the opposing effects of competition at different stages of the invasion 
process: With increasing competition intensity more introduction events are needed until a new 
species can establish, but increasing competition also speeds up the exclusion of the native species 
by an established nonnative competitor. By comparing the ecological and the eco-genetic model, 
we detect and quantify a synergistic feedback between ecological and genetic effects. 
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1. Introduction 



When a new species is introduced to a location where it did not occur before, it begins to interact 
with the resident species, for example as a predator, mutualist, or competitor. These interactions 
are critical in determining the fate of the introduced species and whether and how the community 
changes in response to the introduction. In this study, we investigate the effects of introduced 
species on the native species with which they are competing for resources, such as food or ter- 
ritories. Often there is considerable variation between geographical locations in the impacts of 
an intro duced species o n a particular native competitor, e.g. in container-dwelling mosqu itoes in 
Florida ( jJulianol . Il998l ) or for fish introductions to California (IHerbold and Moyld . Il986l ). Thus 
the question is: Why does the native species suffer from severe impacts or even goes extinct 
in some places, but not in others? While it is evident that the probability for an introduced 
species to become a high-impact invader at a certain location increases with introduc tion rate 
(also known as propaRule pressure, one of the m ost important factors in invasion biology. IDuncanl . 



19971 ; iDrake et all 120051; lLockwood et al.l. 120051) and with th e species' growth rate in the new envi- 



ronment ( Reimanek and Richardson! 1996 ; Wu et al. . 2005), the role of competition has remained 



controversial ( IHerbold and Moyld . Il986l ; IDuncanl . 119971 ; iDavisl . 120031 ) . Some empirical studies sug- 
gest that differences in the intensity of interspecific competition can explain differences in impacts, 
e.g. in the competitive interaction between native bumble bees and introduced honey bees in Cali- 
fornia, where the overlap in the flowers visite d by honey bees and bumble bees was used as a proxy 
for the intensity of interspecific competition ([Thomson! 120061 ) . However, it has been difficult to dis- 
entangle the effect of competition from that of other variables and no consensus has emerged yet. 
Here we contribute a first theoretical building block towards an understanding of how competition 
intensity influences introduced species impacts. 

Considerable impacts of an introduced species on a native species or even extinction of the 
native species can only occur if the introduced species completes two major stages of the invasion 
process: The establishment stage starts with the introduction of one or more founding individuals. 
This new population can then either go extinct, in which case it has to await the next introduction 
event, or increase in size until it is of the same order of magnitude as the competing native 
population. From there, the introduced species can proceed to a second stage in which it becomes 
dominant and may eventually exclude the native species from the community. Thus far, knowledge 
on the effect of competition is limited to single stages of the invasion process. 

At the establishment stage, the intensity of competition with native species appears to have a 
negative effect on the success of introduced species. With ph ylogenet i c rela tedness as a proxy for 
the intensity of competition, this intuitive idea goes back to iDarwinl (118591 chapter 4). Darwin's 
naturalization hypothesis, as it is formulated nowadays, states that introd uced species should be 
less s uccessful at locations where closely related species are already present ( IDuncan and Williams . 
l2002h . Surprisingly, the opposite pattern is somet imes observed i n historical data, for example for 
plants introduced to Hawaii and New Zealand (iDaehlerl . 120011 ; IDuncan and Williams . |2002[ ). A 
possible explanation is that closely related species not only compete with each other but also share 
characteristics that may conf er a high intrinsic ability to survive and grow at the new location 
(IDuncan and Williams! . 120021 ). In the only experimental study on this topic that we know of, 
however, the establishment success of an invader in a microcosm bacterial commun ity increased 
with the average phylogenetic distance to the recipient community (j Jiang et al.l . 120101 ) . In this case, 
it has also been confirmed that more closely related species had a higher overlap in resource use and 
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were therefore competing more intensely. Additional support for the negative effect of competition 
at the esta blishment stage comes from historical bird introductions to Hawaii: iMoulton and Pimm 
fll983f ) and lMoultorj fjl903h observed a negative correlation between establishment success and the 
number of bird species already present, an observation that fits well with the predicti ons of a model 
for the assembly of a competitive Lotka-Volterra community ( iGamarra et all 120051 ) . 

At the second stage of the invasion process, when the nonnative species has already estab- 
lished, the impacts of the introduced species on the native competitors seem to increase with 
increasing intensity of competition. For example, among pairs of congeneric bird species intro- 
duced to Hawaii, pairs in which both species persisted had a significantly higher r elative di f ferenc e 
in beak le ngth than pairs of species in w hich one or both species went extinct ( iMoultonl . Il985l ). 
Contrarily. iRicciardi and Atkinson! (120041 ) and lStrauss et al.l (l2006bl ) found that species less related 
to the native community had higher impacts. However, their measures of impact summarized ef- 
fects on various aspects of the native community, not only those on competitors, such that their 
results do not fully apply to our problem. 

At this second stage, a number of evolutionary and genetic effects can contribute to the im- 
pacts of the introduced species on the native species. If the two species are closely related, there 
is the potential for hybridization and introgression, wh ich can lead to new opportunit ies for evolu- 
tion but also to extinction of rare native populations (jRhymer and Simberloffl . I1996I ). Introduced 
and native competitors also impose new selection regimes on each other which can lead to shifts 
in life hist ories and resou rce use that reduce the intensity of competition and thus facilitate co- 
existence ( Crowded. 119841) . or even a coevolutionary arms race in the exploitation of a limiting 
resource ( iLeger and Espelandl . l2010l ). These evolutionary effects are strongly contingent on the 
species involved. However, there are also genetic effects that are a direct consequence of the reduc- 
tion in native population size due to competition and should therefore be present in most cases: 
smaller populations are subject to inbreeding depression, they can maintain a l ower amount of ge- 



netic diversity, and accumulate d eleterious mutations more rapidly (iFrankhaml . Il995t iLandd . 11995 



Frankham and Kingsolverl . 120041 ) . Additionally, fluctuations in population size or habitat fragmen- 



tation can change the genetic configuration of the native species. So far, there are only a few em- 



pirical studies th at examine such genetic effects of introduced species (e.g. Krueger and Mayl . ll991 



Kim et all 120031 ) . In this study, we focus on the reduction in a native population's genetic diversity 



caused by a population decline after the invasion of a competitor. A reduction in genetic diver- 



sity ca n lower the native species' ability to respo nd to changes in the e nviro nment ( Strauss et al. 



2006al ) and thus lead to a reduced growth rate (lLande and Shannon! Il996l ). This, in turn, leads 



to a further decline in population size, thus closing the feedback loop. Suc h a synergistic feedback 
between ecological and genetic effects can accelerate population extinction ( IRobertl . l201ll ). We will 
call this the eco-genetic effect of the introduced species and quantify how its strength depends on 
the intensity of competition with the native species. 

In summary, a high intensity of competition between introduced and native species appears 
to have contrary effects at the different stages of the invasion process: it makes establishment 
of the introduced species more difficult but also increases the ext inction risk imposed by an al- 
ready established nonnative species on the native species (see also iMacDougall et all 120091 ) . This 
raises the question: What is the overall effect of competition intensity integrated across the entire 
invasion pathway? In this study, we address this question for a scenario in which an unlimited se- 
ries of introduction events would sooner or later lead to the extinction of the native species. Since 
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stochasticity plays an important role in the dynamics of small introduced populations as well as for 
populations on the verge of extinction, we take a stochastic modeling approach and our quantities 
of interest are the expectation and the variance of the time to the extinction of the native species, 
which we will call its persistence time. Very long persistence times in our model can be interpreted 
as signs for indefinite coexistence, since in these cases we would expect ev ents like evolutionary 
divergence of niches to occur before the extinction of the native species (see IStrauss et al. 
and references therein). 

We also study how the relationship between competition and persistence time is modulated by 
the rate at which nonnative individuals are introduced and the nonnative species' intrinsic ability 
to grow and reproduce in the new environment, i.e. we address questions such as: Is the effect of 
competition different for species that are introduced at a high vs. a low rate or does it depend 
on whether the introduced species or the native species has a higher fecundity. First, we consider 
these questions for a purely ecological model, which can be analyzed using theory on birth and 
death processes and a corresponding diffusion approximation. Then we transform this model into 
an eco-genetic model by adding a genetic dimension that allows us to quantify the feedback of 
reduced genetic diversity and adaptability on extinction risk in a variable environment and its 
dependence on competition strength. Lastly, we address a question of interest for invasive species 
management: How low does the introduction rate of a certain species need to be such that the 
native species is expected to persist for a certain threshold time? 

2. Modeling 

2.1. The ecological model 

We represent the population dynamics of the native and the introduced sp ecies as a cont inuous- 



time stochastic model similar to the Moran model in population genetics (IMoranl . 119581 ) . Con- 
sider a community consisting of a fixed number K of individuals, each of which belongs either 
to the native or to the introduced species. The rate at which individuals die is proportional to 
the extent of competit i on ex perienced from conspecifics and members of the other species (as in 



Neuhauser and Pacalal . Il999l ). The strength of interspecific competition relative to intraspecific 
competition is described by the non-negative competition coefficient a. Small values of a represent 
weak interspecific competition and high values intense interspecific competition. In principle, this 
parameter c an be estimated from data on the overlap in resource use between the two species 



( iMayl . [19750 • The individual whose offspring replaces the dead individual is selected by randomly 
drawing one individual from the whole community, including the individual that just died. In this 
draw, native individuals have weight 1 and members of the introduced species weight w. Thus w 
can be understood as the fecundity of the introduced species relative to that of the native species 
in the sense that it is proportional to the per capita number of offspring in a large offspring pool 
from which the new individual is drawn. Note that for the sake of simplicity we assume that the 
competition coefficient a is the same for both species. Thus we are using this parameter to describe 
the symmetric component of the competitive interaction. Nevertheless, an asymmetric situation 
can be generated by setting the fecundity parameter w to a value different from 1. For w > 1 the 
introduced species has a fecundity advantage, for w < 1 a disadvantage. Initially, all K individ- 
uals belong to the native species. From time onwards, single individuals of a nonnative species 
are introduced at rate 7 (on average 7 times per time unit) and start competing with the native 
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species. To bring the community back from K + 1 to K individuals after an introduction event, 
one individual is drawn to die with weights proportional to the competition experienced. 

The population dynamics of the native species can be formulated as a Markov process N = 
{N(t)) t>0 with state space {0, 1, 2, ... , K} which describes the number of native individuals cur- 
rently in the community. Since in this model transitions a re only possible be t ween neighboring 



states, it belongs to the class of birth and death processes ( iKarlin and Taylorl . I1975L p. 131-150). 



The rate A n at which the number of native individuals increases by one is 

c(K — n,n) ■ (K — n) n 
n= K (K-n)-w + n () 

V v ' V v ' 

rate at which probability that a 

members of the introduced r , . . 

species die natl Y e individual 

gives birth 

for n £ {0, . . . , K — 1}, where c(x, y) = x + ay is the competition experienced by an individual 
when the population size of the species it belongs to is x and the size of the competing species is y. 
The constant of proportionality for the death rate is chosen to be 1/ K , such that in the absence of 
the introduced species native individuals die at rate 1. Then there are on average K death events 
per time unit and one time unit can be considered as one generation. For n £ {1, . . . ,K}, the 
rate \x n at which the number of native individuals decreases by one if there are currently n native 
individuals is 

c(n, K — n) ■ n (K — n) ■ w 



K / (K — n) ■ w + n 



rate at which 
native individuals die 



probability that an 
introduced individual 
gives birth 



c(n, K — n + 1) • n 

- ' c(n, K - n + 1) • n + c(K - n + 1, n) • (K - n + 1) ^ 



introduction s v 

ra te probability that a native 

individual dies 



The assumption of a fixed community size is a good approximation for pairs of ecologically 
similar species for which interspecific competition is as strong as intraspecific competition for the 
resource which is limiting population size. This can happen, for example, if there is a fixed number 
of territories or nesting places that can be occupied by one individual from either species. The 
competition for other important resources can be less intense (a < 1) or there can be interspe- 
cific interference (a > 1). The robustness of our model results against violations of the constant 
community size assumption is explored in the supplementary material. In a community of finite 
size, coexistence of native and introduced species is not possible in the long run. We assume that 
our model encompasses the whole range of the native species, such that a reintroduction of native 
individuals from outside is not possible. The introduced species, in contrast, can fail to estab- 
lish and go extinct after an introduction event, but will then be reintroduced at a later time. 
Therefore, the only absorbing state of the model is 0, the state at which the whole community 
consists of introduced individuals and the native species is extinct. Note that in the absence of 
immigration, the native species would be able to persist for an infinite amount of time, since by 
assumption the total number of individuals in the community is constant. In the symmetric case 
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(w — 1), the rare species has an advantage over the more common species for a < 1 because then 
c(n, K — n) = n + a(K — n) < K — n + an = c(K — n,n) if n < K/2. This leads to fluctuations 
of the system around a point (which we will call the coexistence point) in which each species has 
population size K/2. For asymmetric competition (w / 1), whether coexistence is possible on an 
intermediate time scale, depends on the values of a and w. Also the position of the coexistence 
point depends on these parameters. 

Let T n be the random time to extinction of the native species in a realization of the process 
that starts with n native individuals. Let r n and <i\ denote the expected value and the variance 
of T n . We will use the ecological model to compute the expected value tk and the variance a 2 K of 
the time to the extinction of the native species when it is starting with population size K, i.e. in 
the state in which the nonnative species is still absent. 

2.2. The eco-genetic model 

Now we extend the ecological model by including a genetic component. To keep the model 
tractable, we chose the simplest possible genetic scenario: We assume that each native individual 
is haploid and possesses one bi-allelic locus which determines the individual's response to some 
environmental factor, for example whether or not the individual is resistant to a certain parasite. 
At any point in time, one of the two alleles is favored and its carriers have fecundity 1, whereas 
other native individuals have fecundity 1 — s. Thus s is a measure for the strength of selection. 
With probability u an offspring mutates to the respective other allele and at rate e the environment 
and with it the currently favored allele changes. We assume that before introductions start, the 
number of native individuals that carry the favored allele has reached a stationary distribution, i.e. 
it is in mutation-selection equilibrium (see section Appendix B for a derivation of this stationary 



distribution). To be able to compare the results of the eco-genetic model to those of the ecological 
model with the same introduced species fecundity parameter w, we multiplied w in the eco-genetic 
model by the average fecundity w* (see flB.3j) ) of native individuals under the stationary distribu- 
tion. Thus in both models, w can be interpreted as the fecundity of the introduced species relative 
to the average fecundity of native individuals. 

This model can be represented as a Markov process, where the state with n native individuals, 
m of which carry the currently favored allele, is denoted by (n, m). From this state, we can reach 
the states (n+1, m), (n + 1, m + 1), (n, m + 1), (n, m — 1), (n— 1, m — 1), and (n — 1, m) through birth- 
death events, which possibly involve a mutation in the first four cases. The two latter states can 
also be reached through introduction events. The transition rates are defined analogously to those 
in the ecological model, with the death rate proportional to the competition experienced and the 
probability of giving birth proportional to fecundity. The main difference to the ecological model 
is that now the native population is divided into two allelic classes between which individuals can 



switch by mutation. As an example (see Appendix B for all other transition rates), the transition 

(3) 



rate from state (n, m) to state (n, m + 1) is 

c(n, K — n) ■ (n — m) m ■ (1 — u) + (1 — s)(n — m) ■ u 



K m + (1 — s){n — m) + w* ■ w ■ (K — n) 

This is the rate at which one of the n — m native individuals carrying the disfavored allele dies 
multiplied by the probability that it is replaced by a native individual with the favored allele. This 
new individual can either be the offspring of one of the m individuals with a favored allele that did 
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not mutate (probability l — u) or an offspring of one of the n — m parents with the disfavored allele 
that mutated (probability u). The transition from (n, m) to (n, n — m) represents a change of the 
environment and happens at rate e. Our goal here is to compute the expected time to extinction of 
the native species. Here we start in the state (K, m*), where the introduced species is absent and 
m* is the average number of native individuals that carry the favored allele under the stationary 
distribution rounded to the next integer. 



3. Results 



3.1. Ecological effect 

Measured from the start of introductions when the native population has still size K, the 
persistence time of the native species has expectation 




K K i-l , 

1:1: ir 
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1=1 i=i 



0=1 



(5) 



[l + jJiXi 



\ T i+l ~ T i-l, 



iil<i<K 
iii = K 



(6) 



The results (jlj) and (jSJ) are derived by noticing that the time to the extinction of the native species 

when starting in a state n G {1, . . . , K} has the same distribution (denoted =) as the sum of two 
independent random variables: 



T 



Q 1 T 1 ^ O 
— O n + i JV' — &: 




with probability A ^ 
with probability A ^ 



(7) 



for n E {1, . . . , K — 1} and Tr = Sk + Tr~\- Here, S n is a random variable for the time until the 
native population size first changes from state n to a new state N', which in our model is either 
n — 1 or n + 1. S n is exponentially distributed with parameter \ n + fi n for n G {1, . . . , K — 1} and 
with parameter /i K for n = K. Taking t he expectation on both s ides in (J7]) and using T = we 
obtain the following recursion for r n (see lKarlin and Taylorl . Il975l . p. 145-150 for similar problems 
and solutions): 
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if n = 
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if n = K 



(8) 
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which can be solved for tk- Using the law of total variance 



Var(TV) = E [Vax(T N r\N')} + Var(E[T N ,\N'}) = E[<t 2 n ,} + E[t 2 n ,] - E[r N ,} 2 , (9) 
we obtain a similar recursion for the variance 



if n = 

1 I i r 2 \ _| 0n_ ( n 2 — 1_ T 2 ^ 

) 2 A„+/i„ V n+1 1 'n+i; 1 \ n +/i n \ n-1 ' n-lj . . 

\ 2 > (10) 

T^ r n+1 + if 1 < n < K - 1 

which can be solved analogously to (jHJ) once the r n are known (see Appendix C for details of these 



derivations). 

The result given by (j3J) reveals that with increasing strength of competition, the expected time 
to extinction tk decreases until interspecific and intraspecific competition are of similar strength 
(a ~ 1) (Fig. [H A). Here, tk reaches a minimum. If the strength of interspecific competition is 
further increased, tk grows again. The minimizing competition coefficient is below one for low 
introduction rates, and above one for large introduction rates (Fig. [T] B). For low introduction 
rates (7 < 1), the minimum moves towards lower competition coefficients, i.e. weaker competition, 
if fecundities are unequal, no matter whether the introduced species has a higher (w > 1) or a 
lower fecundity parameter (w < 1) than the native species (Fig. [2] A). As expected, tk decreases 
with increasing introduction rate and increasing fecundity advantage of the introduced species. 

These patterns in are paralleled by a corresponding pattern in the variance of the expected 
time to extinction a\. The variance increases with increasing tk and thus also exhibits a minimum. 
To compare the distribution of extinction times to an exponential distribution where the expected 
value equals the standard deviation, we computed the ratio between the standard deviation and 
the expected value of the persistence time in our model. This ratio is close to one for parameter 
combinations that lead to a high persistence time and below one for parameter combinations where 
the extinction of the native species is relatively fast (Fig. [3]). 

The numerical evaluation of (jlj) is practical only for small community sizes. For moderate to 
large community sizes we derived an approximation for the expected time to extinction that is easier 
to compute and gives more insight into the dependence of persistence time on the parameters: 

t k ^K- C — l — f 1 (1 ~ ^ ■ efl«W)-i(i-»0] . e s(t-v) dr] (11) 
Jo (1-U 7 i? V 

where /3 := (a — 1) • K is the rescaled advantage of being common and 5 := (w — 1) • K is the 
rescaled fecundity advantage of the introduced species. 

The result (fTTj) is based on a diffusion approximation of the birth and death process described 
by equations (pQ) and (|2J). The process is rescaled such that the native population size is expressed 
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A 



B 




Fig. 1: The expected time to extinction tk for different values of the introduction rate 7 as 
a function of the competition coefficient a. (B) magnifies the part of (A) around a = 1. The 
minimum of each curve in B is indicated by a solid point. (K = 100, w = 1.) 



as a fraction of the total community and time is sped up by a factor K: 

(N{K-tY 



X 



K 



(x K (t)) 



(^) 



(12) 



f>0 



In the limit as K goes to infinity while and 5 are held constant, Xk converges i n distribution to 
the diffusion process X with infinitesimal generator (see iKarlin and Taylorl . Il98ll p. 195, and the 
supplementary material for a derivation): 



d 



Lf(x) :=-E[f(X(t))\X(0) = x}\ t=Q 



dt 

d 2 d 
x(l — x)-j—^f(x) + ( — /?(! — 2x) (1 — x) x — 7X — Sx(l — x)) —f(x) . 



dx' 



(13) 



The expected time to extinction of the native species whe n its starting frequency is x is a solution 
g(x) of the differential equation (iKarlin and Taylorl . Il98ll p. 193) 



Lg(x) 



1 



(14) 



with boundary conditions g(0) = and | lim ^'(x) | < 00, where \img'(x) denotes the left-sided 

limit at 1 (see supplementary material for details). 

By numerically evalua ting (fTTj) in R ( R Development Core Team . 2010l ) and using a golden 
section search algorithm (jHeathl . |2002| ) we computed r(7,5), the value of /3 that minimizes the 
right hand side of ffTTj) for given values of 7 and 5. After rescaling, we obtained for the competition 
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Fig. 2: The effect of changes in the introduced species fecundity parameter w on the expected 
time to the extinction of the native species and the position of the minimum (A) and on the 
quasi-stationary distribution (d(l), . . . , d(K)) of the Markov process conditional on non-extinction 
of the native species (right column). The curves in the right column correspond to the competition 
coefficients marked with the respective symbol in A. (K = 100, 7 = 0.1.) 
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Fig. 3: The ratio of the standard deviation and the expected value of the time to extinction for 
(A) different introduction rates with w = 1 and (B) for different values of the introduced species 
fecundity with 7 = 0.1. (K = 100.) 



coefficient a which minimizes the expected time to the extinction of the native species. 

a = l+ — ■ (15) 

With increasing introduction rate 7, r(7, 8) increases (Fig. H] A), becoming positive at 7 = 1. The 
absolute value of r(j, 5) increases with the differences in fecundity between the species (Fig. H]B). 

3.2. Eco-genetic effect 

As in (IEJ), we recursively computed T"( n m ), the expected time to the extinction of the native 
species when starting in state (n,m), by decomposing it according to what happens at the first 
jump. Doing this for all states gave rise to a system of (- K+1 )( K+2 ) _ \ li ne ar equations, which we 
solved numerically in R for T^K,m*), where m* is the average initial number of native individuals 
with the favored allele. Although we adjusted the fecundity of the introduced species to match 
the average fecundity of the native species, the expected time to extinction of the native species 
is lower under the eco-genetic model than under the ecological model described in Section 12.11 
(Fig. [5]). Over wide regions in parameter space, T(K,m*) decreases with increasing selection strength 
s acting on the native species (Fig. |5] A) and increases with increasing mutation probability u 
(Fig. |5]B). Fig. |5]C indicates that in the absence of environmental change, the expected time to 
extinction is similar to its counterpart in the ecological model. For a non-zero rate of change, the 
expected time to extinction is reduced, with a particularly strong reduction at small competition 
coefficients. In all cases, the minimizing competition coefficient (indicated by solid points in Fig. 
|5]) is reduced if we take into account the eco-genetic feedback. 
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Fig. 4: The minimizing rescaled competition coefficient r(j, 8) as a function of the introduction 
rate 7 with fixed rescaled fecundity parameter 5 (A) and as a function of 5 with fixed 7 (B). 
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Fig. 5: The expected time to the extinction of the native species under the eco-genetic model as 
a function of the competition coefficient for different selection coefficients (A), different mutation 
probabilities (B), and different rates of environmental change (C). The solid line corresponds to 
the expected time under the ecological model (Eq. 01])). Minima are indicated by solid points. 
(K = 100,7 = 0.1, w = 1.) 
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Fig. 6: The critical introduction rate (on a logarithmic scale) for different fecundity parameters 
w under the ecological and eco-genetic model as a function of the competition strength. In the 
eco-genetic case, w was multiplied by the average fecundity of native individuals under mutation- 
selection balance. (T crit = 2000, e = 0.05, s = 0.25, u = 0.01.) 



3.3. Critical introduction rate 

We set a threshold persistence time r crit and determined t he critical intr oduction rate jcrit, such 
that t k > r crit for all 7 < 'jcrit, using a bisection algorithm (IHeathl . |2002| ) for both the ecological 
and the eco-genetic model. To be able to compare the two models, we adjusted the fecundity 
parameter of the introduced species as above to match the average fecundity of the native species 
under mutation-selection equilibrium in a population of size K. As was the case with the expected 
time to extinction for fixed introduction rate, the critical introduction rate also reaches a minimum 
at an intermediate competition coefficient (Fig. [6]). Not surprisingly, the critical introduction rate 
decreases with increasing fecundity advantage of the introduced species. Taking the ecogenetic 
effect into account, the critical introduction rate is lower than under the purely ecological model. 



4. Simplifications and heuristics 

To intuitively understand why competitors with intermediate interaction coefficients lead to the 
lowest persistence time of the native species, we simplify the state space to three possible states 
(Fig. [7]A): the introduced species is absent and the whole community consists of the native species 
(N), both native and introduced species coexist (N&I), or the native species is extinct and the 
whole community consists of the introduced species (/) . If both species have the same fecundity 
(w = 1) and the introduction rate is small, native and introduced species exclude each other with 
approximately the same probability from the coexistence point and the rate 4>(a) at which this 
happens depends very little on 7. Let ip(a, 7) be the rate at which the nonnative species establishes. 
Then, analogously to ([H]), we can formulate a recursion for the expected time to extinction of the 
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native species: 



T N = 77 T + r Nki (16) 

and 

As solution for the expected time to extinction of the native species when there is currently no 
introduced individuals we obtain 

2 1 

T N 



i[>(a,j) 0(a) ' 

This is essentially twice the sum of the expected sojourn times in the states iV and N&zl. What 
are these times? In the full model, l/i/)(a,j) approximately corresponds to the expected time for 
the introduced species to reach population size K/2 starting from size 0. We will refer to this time 
as the establishment time of the introduced species. The term l/(20(a)) is the expected time for 
one of the two species to go extinct when they are currently coexisting with population size K/2 
each. We will call this the exclusion time. Expressions for establishment time and exclusion time 
as functions of the model transition rates (pQ) and ([2]) were obtained by solving recursions similar 
to equations ([S]) (see supplementary material). 

The expected establishment time is an increasing function of the competition coefficient (dashed 
lines in Fig. 0B and C). The weaker the competition, the higher the advantage of an initially rare 
introduced population and the lower the expected time to reach K/2. The exclusion time (dotted 
lines) on the other hand is a decreasing function of the competition coefficient. The stronger the 
competition, the weaker the force is that drives the system back to the coexistence point, and the 
shorter the time is to the exclusion of one of the two species. Due to these two opposing effects 
the total time (solid lines) to the extinction of the native species, twice the sum of establishment 
time and exclusion time, can exhibit a minimum. The higher the introduction rate, the smaller 
the influence of the competition coefficient is on establishment time, and the flatter the curve of 
l/(20(a)) will be. This is the reason why the position of the minimum is shifted to higher values 
of a as the introduction rate 7 increases (Fig. 0C). At 7 = 1 the boundary where the introduced 
species is absent becomes an entrance boundary for the diffu sion process; this means that the 



process can start at this boundary but can never return to it flKarlin and Taylorl . 1 19811 . p. 235). 
Thus, for higher introduction rates, establishment is no longer a limiting factor. To speed up the 
exclusion of the native species, the minimizing competition coefficient is above one in this region 
of parameter space. In Figs. [2] and H] we observed that in cases where the introduced and the native 
species differ in fecundity, the minimizing competition coefficient differs more from one than in 
the symmetric case. An intuitive explanation for this phenomenon is that in asymmetric cases the 
dynamics are strongly shaped by the differences in fecundity and large changes in the competition 
coefficient are required to affect these dynamics, whereas in the symmetric case small changes in 
the competition coefficient can tip the balance. 

To understand in which states the system spends most of its time, it is useful to examine 
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A B C 




Fig. 7: A simple version of the model with only three states (A) illustrates the antagonistic effects 
of competition strength on establishment time (dashed lines) and exclusion time (dotted line) that 
lead to a minimum in the total time to native species extinction (solid line) which is at a smaller 
value of a in the case of a low introduction rate (B) and at a higher value for a high introduction 
rate (C). (K = 100, w = 1.) 



the quasi- stationary distribution of the Markov process, the limiting distribution of population 
sizes of the native species given that it is not extinct yet. This is computed by eliminating the 
first row and first column, which belong to the absorbing state 0, from the rate matrix of the 
Markov process (see (jA.lft ). The left eigenvector of the remaining matrix associated with the 
eigen value with the largest real part is the quasi-stationary distribution (IDarroch and SenetaL 



19671 ). In the right column of Fig. [2] the quasi-stationary distribution is visualized for different 



parameter combinations. For competition coefficients on the right side of the minimum (indicated 
by squares), the introduced species is absent or has a low population size most of the time. For 
competition coefficients below the minimizing competition coefficient (triangles), the time the 
system spends around the coexistence point contributes most to the expected time to the extinction 
of the native species. In regions of the parameter space where the quasistationary distribution is 
very stable, the ratio of the expected time to extinction and its standard deviation is almost one, 
suggesting that the time to extinction is approximately exponentially distributed. Moving towards 
parameter combinations that lead to a fast extinction of the native species, for example as the 
introduced species' fecundity is increasing (see Fig. [3] B), the standard deviation is decreasing 
relative to the expectation, suggesting that the extinction of the native species becomes more 
deterministic. 

The shift of the minimum in expected extinction time to smaller values of the competition 
coefficient in the model with genetic feedback can also be understood from the quasi-stationary 
distributions in Fig. [2J If competition is weak, the native population is most likely of intermediate 
size, and therefore a substantially lower amount of genetic diversity can be maintained within the 
population compared to a population that makes up the whole community. For high competition 
coefficients, the native population is most likely near its carrying capacity, if it is still present, which 
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we condition on, and has therefore almost its full adaptability. Thus the additional extinction risk, 
or the reduction in expected time to extinction, would be larger for small competition coefficients 
than for large competition coefficients. This effect can be seen in Fig. |5l as the curves for the 
eco-genetic model diverge more from the curve belonging to the ecological model at lower than 
they do at higher competition coefficients. 



5. Discussion 



Our theoretical results indicate that the introduction of nonnative competitors raises the ex- 
tinction risk of native species, both directly and indirectly, via a reduction in genetic diversity. 
The expected impact does not generally increase or decrease with competition intensity as one 
might expect, but there is an intermediate competition coefficient for which the expected time to 
extinction of the native species is minimized. This is the result of the opposing effects of compe- 
tition strength on the establishment step of the invasion process and on the impact of an already 
established species. Introduced species that do not compete intensively with species from their new 
range can readily establish, but their ecological impacts are weak and it will take a long time for 
them to drive one of the native competitors to extinction. On the other hand, an introduced species 
that is competing very intensely with one of the native species has high potential ecological impacts 
once it has established. However, such a species may need a lot of introduction attempts before it 
can establish, because the native competitor can efficiently exclude it from the community. Based 
on our results we expect competitors with intermediate interaction strength to also have the lowest 
critical introduction rate. Thus if one would set a management target to preserve an endangered 
native species for a certain time period, then the greatest introduction prevention efforts would be 
necessary for nonnative species that would have an intermediate intensity of competition with the 
native species. 

With the help of our eco-genetic model, we quantified the feedback between ecological and 
genetic effects of the introduced species on the native competitor. This feedback is synergistic in 
the sense that ecological and genetic effects enhance each other: A reduction in population size 
causes a reduction in genetic diversity and this reduction in diversity can lead to further popu- 
lation decline in a changing environment. We found that this eco-genetic feedback is particularly 
strong for small intensities of competition between the introduced and the native species. This 
is because for high competition coefficients, the native species spends most of the time until its 
extinction in states with a high population size, whereas for intermediate and low competition 
coefficients the native species spends most of the time coexisting with the introduced species at 
an intermediate population size with a corresponding reduction in genetic diversity. Because the 
eco-genetic effect makes competitors with a relatively small competition intensity more dangerous 
for the native species, the minimizing competition coefficient is smaller compared to the ecological 
model. Similarly, also the critical introduction rate is lower than in the purely ecological scenario. 
This highlights the importance of including eco-genetic feedbacks into risk assessment models. If 
we only bring the introduction rate down to the level required under the ecological model and 
there is an eco-genetic feedback, we will fail to keep the expected impacts below the prescribed 
threshold. Eco-genetic effects are also a possible explanation for why there is relatively little evi- 
dence in the invasi o n bio logy literature for native species extinction due directly to an introduced 
competitor (IDavisl . 120031 ): Most endangered species are not threatened by a single stressor, but 
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by co mbinations of them, for example habitat change and invasive species ( iGurevitch and Padillal . 



2004 ). The introduction of a competitor can weaken a native species' ability to respond to other 



stressors. There may have been many extinction events which were attributed to other factors and 
in which a significant contribution from an introduced competitor went unrecognized. Thus, when 
making predictions on native species population dynamics it can be important to c onsider the pos- 
sibly synergistic interaction of species invasions with other drivers of global change (jDidham et all 



2003). 



Previous theoretical studies on the impacts of introduced species have buil t models designe d 



to understand these impacts in specific systems (IByers and Goldwasserl . l2001t iThomsonl . 120051 ). 



Here we contribute a building block towards the development of a theory that predicts impact 
from parameters of the introduced species, the native community and the introduction process. In 
this first stage of theory development our focus has been on simple models that are analytically 
tractable and give us insight into general phenomena. Of course, these models could be extended in 
many ways to incorporate more biological realism or to adjust them to specific biological systems, 
e.g. by including age or stage struc ture, which might in some cases influence the outcome of invasion 
and extinction dynamics (see e.g. lLande and Orzackl . Il988l ). Nevertheless, we conjecture that our 
main findings are robust to several kinds of changes in the model formulation. As an example, if 
one would assume that competition affects fecundity instead of viability, the transition rates of 
the Markov process would look slightly different. For the diffusion approximation, however, only 
terms that are linear in the parameters play a role such that it would remain unchanged. Even a 
model in which the assumption of a fixed community size has been relaxed converges to the same 
diffusion approximation (see supplementary material). Since the expected persistence time under 
this diffusion approximation exhibits a minimum at an intermediate intensity of competition we 
know that this phenomenon must also be present in these and a large number of other modified 
models. 

In this study, we provide a model for the expected effect of an introduced species on one 
native competitor. Of course, native communities may consist of multiple competitors, as well as 
predators, mutualists, and parasites. To predict the impact of an introduced species on a whole 
community, our model could be combined with models for the other components of the community 
and interactions be tween them . Such detailed models have not been analyzed thus far. However, 
Fig. 2 in a study by lCasd (1l990l ) shows that the probability that an introduced species can establish 
in a multi-species competitive community and replaces at lea s t one native species is maximized 
for intermediate mean competition coefficients. Although ICase fll990l ) does not address this point, 
this is one hint that our finding that intermediate levels of competition are most dangerous scales 
up to more complex communities. 

In our model, single individuals of the nonnative species were introd uced into the new h abitat. 
What happens if multiple nonnative individuals are released at once? iDrake et al.l (120051 ) found 



that the product of introduction frequency and introduction size was a good predictor for the 
persistence of introduced Daphnia populations and that adding introduction frequency and size 
as single factors did not lead to significant improvements. However, if the nonnative population is 
subject to an Allee effec t, i.e. positive density- dependence of population growth at low densities 
(ICourchamp et allll999l ). the expected establishment success can strongly differ between a scenario 
with freque nt introductions of one or a few individuals and one with rare introductions of many 
individuals (IDrury et all 120071 ). Disentangling the effects of propagule size and propagule frequency 



17 



for such important scenarios is a promising field of future research. 
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Appendix A. Rate matrix of the ecological model 

Ecological model 

The rate matrix of the ecological model is 
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/-A A 

Mi -(Ai+^i) Ai 
\i2 -(A 2 +A*2) A 2 
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A* if -^k ) 



(A.l) 



where the Ai and /ij are given by equations jT]) and ([2]). 

Appendix B. Full specification of the eco-genetic model 

Let (n, m) be the state with n native individuals, m of which carry the currently favored allele. 
Then the transition rates for 1 < n < K and < m < n are 
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(B.l) 

where 

obtain w*, the average fecundity of the native species under the eco-genetic model in the absence 
of the introduced species {n = K), we computed the stationary distribution d = (<i(0), . . . , d(K)) 
of the birth and death process that describes the number of favored alleles in a native population 
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of size K and has transition rates 

'm + 1: (K - m) ■ Mi-u)+(i-s)(K- m )u for ^ 

\ / m+(l— s)(iv — m) 

m — 1 : 



m — y < 



m+(l— s)(-ff-m) 
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— — — - '- tor m > 
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K — m : e for < m < K 



(B.2) 



and averaged 



K 



w 



(B-3) 



Appendix C. Recursive solution for expectation and variance of the time to extinc- 
tion in the ecological model 



The middle equation in the recursion ([8]) can be rewritten as 

A*n ( r n — r n-l) — 1 + ( r n+l — T n) ■ 

Define z n := r„ — r n _i, such that 



Solving this recursion for z with 2;^ = — r^_i = — gives: 



K i-1 

E^n^ 



A*i A*j 



(CI) 



(C.2) 



(C.3) 



m K ^ i—1 



E* = EE^- 

1=1 1=1 i=l ^ l j=l ^ 



With m = K we obtain the expected persistence time (jJJ). 
Similarly, the middle equation in ( flu]) can be written as 

c — 1 ^ n c 
/in 



(C.4) 



(C.5) 



where = 0^ — o\_ x and r]„ is given by (jSJ). As above, we can solve this recursion for £ with 

6r = 



K i-1 



(C.6) 
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and finally obtain (jHJ) by summing 

K 
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Supplementary material 

Diffusion approximation to the ecological model 

We are seeking an approximation of equation (J4j) for large community sizes K. The strategy 
used here is to first rescale the Markov process defined by the transition rates (pQ) and ([2]) and 
approximate it by a diffusion process. Then we compute the expected time to native extinction 
in the diffusion process and finally prove that the exact recursive solution given by equation (j4]) 
converges to the expected time under the diffusion process as K goes to infinity. 

Convergence in distribution X^ =>■ X 

Theorem 1. As the community size K goes to infinity, the rescaled Markov process 

x K = (x K (t)) t > = (^§^) t>o - (s- 8 ) 

with a = 1 + -4 and w = 1 + ~, where j3 and 5 are constants, converges in distribution to a 
diffusion process X with infinitesimal generator 

Lf{x) = \ Kx) o^ f{x) + a{x) Tx f{x) ' (S - 9) 

where 

a(x) = — (3(1 — 2x){l — x)x — 5x(l — x) — jx (S.10) 

is the infinitesimal mean of the diffusion process and 

b{x) = 2x{l-x) (S.ll) 

is the infinitesimal variance. 

To prove the theorem we need a few lemmata. 
Lemma 2. For x G {0, j^, . . . , 1} the generator of the Markov process X K , defined as 

L K f(x) := jV[f(X K (t))\X K (0) = x] \ t=Q , (S.12) 

where f is a bounded twice continuously differentiable function, converges to the generator of the 
diffusion process given by equation IIS. 9\) : 

lim max \L K f(x) - Lf(x)\ = . (S.13) 



Proof. 
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(S.14) 
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f ' {x) + l f " {x) --R + <D \Ki 



Xrx + 



fj, Kx (S.15) 



f{x) ■ (X Kx - vk x ) + \f"{x) ■ Xk * + ^ k * + O (± ) • (X Kx + hk* 



(S.16) 



Substituting the scaled parameters into equations §B) and and using Oj(l) to denote terms 
which fulfill lim Oj(l) = uniformly in x G [0, 1], we obtain 



K— too 



Xrx — fJ-Kx —K 



_ 7 . 
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l + o 6 (l) 



Thus 



L K f{x) = x(l-x) 



dx 2 



(S.18) 



(S.19) 



f(x)+ (-p(l-2x)(l-x)x-ix-6x(l-x))—f(x) + o a (l) . (S.20) 

dx 



All expressions are bounded uniformly in {^,^, ■ ■ ■ ,%} and the error in equation ( IS. 131) thus 
converges to zero uniformly in { , , . . . , ^ } as A' goes to infinity. □ 

Lemma 3. The sequence of Markov processes (Xk) KgN is tight. 

Proof. This is proven by using the compactness o f the state space and the basic criterion for 
tightness and the Aldous condition from p. 34-35 in lJoffe and Metivierl ( 119861 ). □ 



Lemma 4. The martingale problem for the generator L has at most one solution. 

Proof. The problem can be written as a stochastic differential equation with infinitesimal parame- 
ters extended beyond the interval [0, 1] with a(x) = a(0) and b(x) = b(0) for x < and a(x) = a(l) 
and b(x) = 6(1) for x > 1 such that the process is now defined on K while leaving the behavior 
inside t he interval [0, 1] unchanged. Then we can apply the Yamada-Watanabe theorem (Theorem 
26.10 in lKlenkdl2008l ) and conclude that the stochastic di fferential equ ation has a unique strong so- 
lution. This implies uniqueness in law (Theorem 26.18 in iKlenkd 120081 1. which in turn is equivalent 
to the statement that there exists at most one solution to the corresponding martingale problem 
dRogers and Williams! . I2OO0L p. 159). □ 
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Proof of The oremUl The conve r gence in distribution follows from lemmata EJ [31 and H] and Theo- 
rem 4.8.10 in HthierandKurt3 (J2005|). □ 



Expected time to extinction under the diffusion process 

Now we compute the expected time to the extinction of the native species under the diffusion 
process. Define 



ai(X) := inf{t > : X(t) < 1} 
and let E x denote the expectation for a process starting in x. 
Theorem 5. 

(l-0-7ln(l-C)+^ ' e 



(S.21) 



E x [a (X)}= I i 
h 



77(1-77) 



dr) d£ 



(S.22) 



for x E [0, 1]. 



Proof. An important tool for the analysis of diffusion processes is the scale function S(x). The 
scale function transforms the state space such that the diffusion pr ocess becomes a marting ale. 
Thus the scale function fulfills the differential equation LS(x) = ( Karlin and Taylor] . 19811, p. 
196). With s(x) := j-S(x) the differential equation becomes 



-b(x) — s(x) + a(x)s(x) = . 
2 ax 



Consequently, 



S(X) 



:, J() 



x _ 2a(y) 

b{y) 



dy 



Jo ' -[-P(.l-2y)-^-S]dy = J3x(l-x)-~fln(l-x)+8x 



(S.23) 



(S.24) 



where the lower limit o f the integral can be chosen arbitrarily and is here for convenience 

(IKarlin and Taylor! . Il98ll . p. 194). With this we can write the infinitesimal generator as (IKarlin and Taylor . 
198l|, p. 195): 



s(x) 



(S.25) 



Let g(x) be the solution of the differential equation (IKarlin and Taylorl . Il98ll p. 193) 

Lg(x) = -1 (S.26) 

with boundary conditions g(0) = and | lim(yf'(x)| < 00. Using the differential operator from 

x yi 

(IS.25p . equation (1S.26[) is equivalent to 



e /3f (1-O-7 Mi-0+<5£ (j i _ 
V ^0.5 



€ g-^ r ?( 1 - f ?)+7l n ( 1 -»?)- 5 »? 



77(1 - 7]) 



drjj d^ + C 2 



(S.27) 
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The boundary conditions are fulfilled with C2 = and 



1 e — /3r^(l— 77)+t ln(l— 17)— <5?7 



Ci= L .a-,) (s - 28) 



With this 



x _ _ _ /•! e - / 37?(l-j?)+7ln(l-j?)-5r? 



^(x)=/ e «d-0-7Mi-0+« / r (s.29) 

Jo ^(1-^) 

Since g(a;) < g(l) < 00 for all x G [0, 1], g is a bounded, twice continuously differentiable function. 
Then, since (Jf (i)) t > solves the martingale problem for the generator L 

g(X(t))- [ Lg(X(s))ds = g(X(t)) + t (S.30) 
Jo 



with t e [0, 00) is a martingale. Application of the optional stopping theorem (see lEthier and Kurtzl . 



20051 . p. 421, for a similar application) with r > gives 

E x [g (X(a (X) A r))] + [a (X) A r] = g(x) . (S.31) 
Letting r go to infinity, we obtain by the dominated and monotone convergence theorems 

g{x) = MX)] . (S.32) 

□ 

Convergence of extinction times 

Because the first time to reach the boundary is not a continuous functional of the process, we 
cannot conclude from the weak convergence of the rescaled birth and death process to the diffusion 
process that also the expected time to extinction converges to that under the diffusion process. 
Some more work is required to prove 

Theorem 6. 

lim ExMXaO] = E 1 [a (X)] . (S.33) 

K— >oo 



We start with proving the following lemmata: 
Lemma 7. For all j E {1, . . . , K — 1} 



(S.34) 
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Proof. Substituting the scaled parameters given in theorem CD into equations (DO) and (J2J), we find 
that 



A, 



and 



(K+%-j)-(K-j)-j 
K- [K+i-(K-j)] 



[K+%.(K-j)].j-{K-j)(l + ±) 



+ 7 



(S.35) 



[K + l + ^.(K-j + l)] -j 



(S.36) 



[K + l + &.(K-j + l)]-j+[K + l + &.j].(K-j + l) 
From equations IS. 351 and IS. 361 it follows that 







+ 7 - 


[1 + ^ + 




j 

K 


+ *)] [1 + 1(1- 


*)] 


Xj 1 + 


j 

K K 


K-j 


.(l + i) S 


i- ^ K 


(1 


_ ± + ±) £ 

K ^ K) K 


(i + 


j \ 
K K ) 



[l + #(l-i)](l + D , 7 

After factorization one obtains ( IS . 34[) . 
Lemma 8. 

2 . (1 ~ ikY 



i+#i v U 



i_1 A 

n-< 



(i-*) 7 



• e°« < e^'+l 5 ' • e°« =: E 



(S.37) 



(S.38) 



□ 



(S.39) 



for all l,i G {1, . . . , i^} wrai/i / < z. 
Proof. Using lemma [7] we can rewrite 

n- 



as 



exp 



i-l 



l + #(l-i) 



1 — K K 



1 + T7 + 



7 



if K-j 



exp 



i-l 



E ln ( 1+ I^ 



3=1 



K 



In 1 + — 1-^7 - In 1 + — + 



K 



7 



K K-j 



(S.40) 



. (S.41) 
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Since ln(l + x) = x + 0(x 2 ) for x close to zero and — ln(l + x) < —x + x 2 for all x G (—0.5, oo) 

i-l 



3=1 



J 


H) 


_S_ , 1 


KK K 




K K - j (K - j f 



and 



°4 



(S.42) 



i-l 



i=z v J; i=i J 



With this 



i-l 



\p J p 


H) 


5 7 


KK K 




K K-j 



0(1) 



(S.43) 



(S.44) 



<2|/3| + |5|-7ln n?—L\+0(l) 

and inequality (IS.39[) follows from this. 

Lemma 9. There exists a function r(ip) such that 

lim EVM-Xx)] <KV0 

K— >oo 

and 



(S.45) 
□ 

(S.46) 



lim r(^) = . 



(S.47) 



Proof. 



K K 



^E^E*< + 

i=i i=i 



Z[Kil)~\ 

K 



(S.48) 



with m = \Kip~] — 1 and where \x] denotes the smallest integer larger than or equal to x. 7$ is 
given by equation (1C.4j) and Zi by equation (1C.3j) . We first consider the first summand: 



m K ^ i—1 



;=1 



(S.49) 



1=1 i=l 



i i—1 



K—m m i—1 



K-l 



m i—1 



m K-l 



_, / lib b b _L» J \ lib _, lib b _L ■» 11 _L lib b _L» _, I I b ±\. _L» \ 

f E^n^E^En^ e ^En^^En^ 

A \ i=l ^ 1=1 j=l ™ i=m+l to 1=1 j=l to i=K _ m+1 to l=1 j=l fJ-K l=1 j=l 



A 



B 



C 



D 
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Note that 



m > k- 



With this 



K \ K 



[! + *(!-*)] ( 1 + i)_ >K . i 



K V K 







1 if 





. (S.51) 



E^ 



K ~ K' ^ iii- HK ^ 1 - 4 - //X 1 - f - if 1 - -0 ' 

i=l i=l K 



Em E ib 
< < 



(S.52) 



and 



B <m -E K < 2m • E K" _ 4m • E f K \ < 4ib ■ E n f 1 

^ i(K-i) ~ HK ^ 7^4 ~ ///\ 11 I 2m" 

2 



K ~ HK 



i=m+l 



H 



In ( - ) , (S.53) 



and 



_E • m 



C / '. v -> 1 v • \ i 

^ £ ™ ,.ik» ^i) Hf £ HK • (A' - m) • (1 - ff 



E i 



m+1 



7-1 



(S.54) 



Let /(x) = (1 - x) 7_1 . Then f'(x) = (1 - 7 )(1 - x)^ 2 . f'(x) < if 7 > 1 and f'(x) > if 7 < 1 
for all x E (0,1). If f'(x) > 



K-l , . 

E f {lc 

i=K—m x ' 

is a lower sum of the integral 
K ■ [ f(x)dx 

and if f'(x) < 



(S.55) 



(S.56) 



K 



E / 



(S.57) 



i=K-m+l 

is a lower sum. Thus, in both cases 



K-l 



i=K-rn+l 



7-1 



K /m\T 



E <*/ m <i-*r'*=^)'. 



(S.58) 
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Consequently 



C E-m 1 /my E ■ 

K ~ H ■ (K — m) ■ {l — j±y 7 \KJ ~ H • 7 • (1 - ^) 7+1 ' 1 J 

Since 

(# + 1 + i) ■ K 7 
u x = 7 • V ^ ^ > 1 , S.60 

^<^<^. (S.61) 
K ~ K ■ fi K ~ 7 1 J 

Turning to the second summand in equation (IS.48I) and using equation (10.3() 

k ^ i-l 

z \Kif\ 

< £ v 1 , g < £ in(^-r^i) Jg r , 

- H J-£ i-(K-i) HK-K ~ H \K1>] 7 • K l^o ■ 1 J 



In conclusion 



lim EJa (X K )} < — ( — ^— + #ln f — ) + - 

^ <n ^ ^2^/ 7-(l-V) 7+ V 7 

To see that lim r(ib) = note that 



1 \ , ■ £> ( ^ J 



lim ^ ■ In [ — = lim — — = lim 4- = lim -0 = 0. (S.65) 



Proof of Theorem 



□ 



lim |Ei [<T {X K )) - Ex [<t (X)]| (S.66) 

A— >oo 

lim |Ei - Ei [<ri{X K )\ + Ei [oi(Xjr)] - E 1 [^(X)] + E x [a t (X)] - E 1 [a (X)}\ (S.67) 

A • X 



< lim Ex [a (X K )\ - E x [<t,(X*)] + |E X [^(X^)] - E x [<r,(X)]| + |E X - E x [a (X)]| 



(S.68) 
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<r(l)+ lim \E 1 [(x l (X K )]-E 1 [a l (X)]\ + \E 1 [c7 l (X)]-E 1 [a (X)]\. 



K^oo 



(S.69) 



provided that the limit in the second summand exists. Following Kurtzl ( 1981 ) p. 13-14 we argue 
that for a fixed path x, o~i(x) is decreasing in / and can therefore be discontinuous at only a 
countable number of values of /. Since the diffusion process is a rescaled Brownian motion, there 
is no point in the interval (0,1) that is exceptional compared to the other points and the points 
of discontinuity are placed according to some continuous probability distribution on the interval 
(P (discontinuity at I) = V/). For any discontinuity point thus the collection of paths for which 
the discontinuity is placed exactly at I is a null set. Since the union of countably many null sets 
is a null set itself, the probability that a random path has a discontinuity at I is zero. Therefore, 
the weak convergence of Xk to X implies the co nvergence of <ti(Xk) to &i{X) according to the 
Continuous Mapping Theorem (Theor em 13.25 inlKlen ke 2008) such that the second limit exists 
and is equal to 0. See Lemma 3.3 in IChigansky and Klebanerl ( 120121 ) for a more rigorous proof 
requiring a positive and increasing scale function S(x), which we can achieve by choosing 



S(x) 



3 +/3y(l-2/)-7 ln(l-j/)+«5i/ fly 



(S.70) 



and a positive quadratic variation [X, X] t for t > 0. 

After the middle summand vanished, we take the limit Z — >- 0: 

lim |Ei [a (X K )} - E x [a (X)]\ < limr(Z) + lim {E^X)] - E 1 [a (X)]\ = + = 0. (S.71) 

The first limit follows from Lemma |9j the second limit from the boundedness of expected times 
under the diffusion process (Theorem [5]) and the application of dominated convergence. □ 

Consequently, the expected time to extinction under the diffusion process is a valid approxi- 
mation for the expected time under the rescaled birth and death process. Figure IST81 visualizes the 
quality of the approximation for different community sizes. 

Derivation of establishment time and exclusion time 

The recursion for the expected time f n to reach K/2 from a population size n above K/2 is: 







\n~\-f-ln An 



T n +1 + 



if 71= f 
K 



a„+^-i iff + l<n<K-l . 
if n = K 



(S.72) 



This can be solved in the same way as recursion (JHJ). For the expected establishment time if the 
native species starts from its carrying capacity, we obtain 



K K 



i-1 



l=K/2+l i=l ^ % j=l ^ d 



(S.73) 
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Fig. S.8: The exact solution for the expected time to the extinction of the native species (jlj) 
approaches the diffusion approximation (1S.22j) as the community size increases. (7 = 0.1). 



We compute the exclusion time Tkii (the time to reach either or K from K/2) under the 
assumption of symmetry and set 7 = in model equation [2j The appropriate recursion is: 



if n = 

t— j; — + T^f n+ i + T^fn-i iil<n<K-l . 

if n = K 



(S.74) 



Equations fIC.ip and (10. 2[) remain valid but now we have Z\ = f\ — r = f\ and zk = t~k ~ t k-i 
—t k _i, such that 



K-l > K-l i-1 , 

-^n-+E 1 n- and 

j=m J i=m 3=m J 



(S.75) 



K-l 



K-l „ i-1 



A, 



n=z 1 = -tk-i n — + yi - n 

j= i Pi i= i j=1 H 



(S.76) 
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Because of symmetry 



£—ii=\ Ui 1 1 7 = 1 u, 



£—ii=l Ui 1 1.7=1 u,- 



i + n^ -1 ^ 



and 



(S.77) 



Under the diffusion approximation, expressions for establishment time g and exclusion time g 
can be found by using appropriate boundary conditions in equation (1S.27|) . For the establishment 
time in the symmetric case, we need to use g(l/2) = and | lim^'(x)| < oo, such that 

9(1) = f tt^t- f (1 " T7)7 " 1 • e^-^drj d£ . (S.79) 

For the exclusion time, we need to use the boundary conditions g(l) = g(0) = and we set 7 = 0, 
which leads to 

/■1/2 P -Mi-v) rv , 

~ 9i °- 5) = J V (l - i) I 6 ^ dV ■ (S - 80) 

The assumption of a fixed community size 

If two competing species differ in their specialization on a range of resources then we would 
expect the habitat to support a higher total number of individuals when both species are present 
than if either of the species is alone and the other one absent. Let's assume that the two species are 
competing with competition coefficient a and have the same intrinsic growth rates r and carrying 
capacities K in the absence of the other species. Then the dynamics of the system can be described 
by two coupled ordinary differential equations (a case of the classical Lotka-Volterra competition 
model) 

dNi ( Ni + aNA „ .„ . 

r • Ni ■ 1 °- for i = l,j = 2 or i = 2, j = 1 (S.81) 



dt V K 

Because the situation is entirely symmetric, the equilibrium population sizes of the two species 
should be identical: = iV| = N*. We can compute it by solving 

/ N* + aN*\ 

r • iV* • ( 1 - ™ J = . (S.82) 



Thus 

2-N* = 

1 + a 



2K 

2-N* = - (S.83) 
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is the total size of the community at equilibrium. In the derivation of the diffusion approximation 
( llljl . we assumed that as K goes to infinity a approaches 1, such that 



(a- 1) -K^p. (S.84) 

If this holds, then also 
IK 

(a - 1) • — > . (S.85) 

1 + a 

Thus, whether or not the total size of the community is held constant or follows the Lotka-Volterra 
competition model, we obtain the same diffusion approximation as long as the advantage of the 
rare species is sufficiently weak such that the scaling equation (1S.84|) is fulfilled. 
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